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Cancer metastasis is an important criterion to evaluate tumor malignancy. Matrix metalloproteinases 
(MMPs) play a crucial role in cancer proliferation and migration by virtue of their proteolytic functions in 
angiogenesis and extracelluar matrix (ECM) degradation, making them potential targets of anti-metastaic 
therapeutics. Recently we showed with both in vivo and in vitro experiments that metallofullerenol Gd@ 
C82(OH)22 can effectively inhibit MMP-2 and MMP-9 with high antitumoral efficacy. Furthermore, our in 
silico study revealed that Gd@C82(OH)22 could indirectly inhibit the proteolysis of MMP-9 via allosteric 
modulation exclusively at the ligand specificity SI' loop. Here, we expand our study toward another 
gelatinase, MMP-2, using molecular dynamics simulations. Despite the high structural similarity with 
64.3% sequence identity, their responses to Gd@C82(OH)22 were quite different. Toward MMP-2, Gd@ 
C82(OH)22 could block either the Zn^'''-catalylitic site directly or the 81' loop indirectly. Surface 
electrostatics uniquely determines the initial adsorption of Gd@C82(OH)22 on MMP-2, and then its further 
location of the most favorable binding site(s). These findings not only illustrated how the inhibitory 
mechanism of Gd@C82(OH)22 is distinguished between the two gelatinase MMPs with atomic details, but 
also shed light on the de novo design of anti-metastatic nanotherapeutics with enhanced target specificity. 

Over the past decade, nanomedicine has emerged as an alternative and novel approach for overcoming 
some of the difficulties experienced by conventional medical approaches in delivering the benefit of 
small molecule-based medicines' \ There are many applications of nanomedicine, such as drug delivery, 
diagnostics and therapy. Nanomaterials offer a wide range of applications because of their unmatched flexibility 
and diversity, endowed by the wide selection of materials including polymers, liposomes, inorganic materials, and 
carbon alIotropes'''\ Another unique property of these nanomaterials, the intrinsically large surface-to-volume 
ratio, allows multiple biochemical and functional groups or "blocks" to be attached onto their surfaces'", making 
them effective templates for de novo nanomedicine design to overcome the limitations often exist in conventional 
drugs, such as low target selectivity and poor release control. Many different platforms have been devised and 
specified with single or multiple functions aiming at applications for target selectivity, payload optimization, and/or 
release control. Roughly speaking, nanomedicines can be categorized into diaganostics and therapeutics depending 
on their general application purposes. Based on their role in pharmacokinetics, therapeutic nanomedicines can be 
further categorized into nanocarriers for conventional drugs and nanodrugs for direct treatment of the target 
diseases. There have been considerable efforts in the research of nanocarriers, such as the FDA approved Abraxane 
and Doxil, but relatively fewer studies on direct nanodrugs (hyperthermia' and photodynamic therapy^ might be 
loosely counted into this category). Despite these efforts (and promising achievements), there is still a lack of 
detailed understanding of molecular interactions between nanomaterials and biomolecules, which is crucial for 
developing better and smarter nanomedicines, as well as avoiding undesirable side effects"". 

Meanwhile, gadolinium ion (Gd""^) has been extensively used as contrast agents for enhanced magnetic 
resonance imaging (MRI) due to its high proton relaxivity'^ ". Recently, fullerene Cs2 was shown to serve as 
an effective host to safely secure the toxic Gd"** with even better contrasting power than conventional agents after 
water-soluble hydroxylation (i.e., Gd@C82(OH)22)'''''^ More importantly, Gd@Cs2(OH)22 was then found to 
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have strong anti-tumoral efficacy, effectively inhibiting tumor growth 
and metastasis in human hepatoma and breast cancers xenographted 
in mice". Even at an ultra-low dose, Gd@C82(OH)22 shows compar- 
able anti-tumoral efficacy with conventional antineoplastic agents 
(e.g., cyclophosphamide and cisplatin) without noticeable side effects. 
Gd@C82(OH)22 also circumvents the acquired resistance to cisplatin 
by reactivating the impaired endocytosis of human prostate tumor 
cells'^ 

Subsequent study showed that Gd@C82(OH)22 simultaneously tar- 
gets more than 10 angiogenic factors by down- regulating mRNA and 
protein expression levels, resulting in a dramatic decrease of tumor 
microvessel density (MVD), thus lowering blood supply to the tumor 
tissue". Very recently, we reported that Gd@Cs2(OH)22 effectively 
suppresses volume growth of human pancreatic tumor xenografted 
in mice by >50% compared to the saline control group'". The envir- 
onmental scanning electron microscope (ESEM) shows clear and 
smooth surface in the tumor tissue treated with Gd@Cg2{OH)22, 
which is again highly contrasted with the saline control one where 
the typical rough surface is heavily and densely packed with blood 
vessels. Anti-CD31 antibody stains also showed evidence of fewer 
MVD and lower CD31 expression in the Gd@C82(OH)22-treated 
group, indicating restriction of angiogenesis. More specifically, we 
also identified that Gd@Cs2(OH)22 effectively inhibits the gelatinases, 
matrix metalloproteinase-2 (MMP-2) and MMP-9, in both mRNA 
expression and enzyme activity levels"*. Under normal physiological 
condition, MMPs are tightly regulated, as they govern the proteolysis 
in various contexts such as in tissue remodeling and wound healing''^ 
Once recruited in the pathological environment such as malignant 
tumors, however, MMPs are involved in new blood vessel formation 
and ECM degradation^"'^', modifying angiogentic ECM^^'^' and cell 
surface molecules^'', as well as mediating the release of ECM-bound 
pro-angiogenic factors like VEFG and FGF-2^^'^''. Therefore, MMPs, 
specifically their Zn^^-coordinated active sites, have been designated 
as a potential target for anti-metastatic therapeutics^'. However, sim- 
ilar to kinases^** the structural similarity in the catalytic sites among 
all MMPs (at least 28 of them in total) often limits the target selectivity 
even with favorable binding affinity, resulting in undesired side 
effects". 

Furthermore, our in silico study with molecular dynamics simula- 
tions revealed that Gd@C82(OH)22 could favorably interact with 
MMP-9 at the ligand specificity SI' loop, possibly regulating the 
ligand binding onto the catalytic site remotely'". Interestingly, the 
SI ' loop is known as the most characteristic region from the catalytic 
domain in terms of length, sequence similarity, and ligand recog- 
nition among the available MMPs". Thus, our finding, in contrast 
with conventional MMP inhibitors (MMPi's), proposed an alterna- 
tive inhibitory pathway that Gd@Cs2(OH)22 allosterically modulates 
the native ligand binding on MMP-9 via an indirect but very specific 
exocite interaction. From a different perspective, this also demon- 
strates the possibility of applying the carbon-based nanoparticle as a 
direct nanodrug through the specific interactions with the metastatic 
protein beyond the general role of nanomedicine as a nanocarrier for 
conventional drugs'*^. 

In this work, we extend our investigation to the interactions between 
Gd@C82(OH)22 and the gelatinase MMP-2 using atomistic molecular 
dynamics simulation. In our previous in vitro studies, activities of both 
MMP-2 and MMP-9 were inhibited by Gd(a)C82(OH)22, even though 
MMP-9 seems to be more dramatically affected (Fig. SI)"*. Initially, we 
presumed that the binding dynamics of Gd@C82(OH)22 would be 
similar for both gelatinases, since MMP-2 shares a highly conservative 
catalytic domain structure with MMP-9 (i.e. RMSD = 0.97A) (Fig. 1). 
However, our molecular dynamics simulations showed that MMP-2 
exhibit quite different binding modes with Gd@C82(OH)22 from 
MMP-9, thus implying different inhibitory mechanisms between 
them. The underlying deviations in their molecular binding mechan- 
isms are also illustrated in detail. 



Results and discussion 

The molecular system was set up with the catalytic domain of MMP- 
2 (PDB code: IQIB)" surrounded by multiple nanoparticles of metal- 
lofuUerenol Gd@C82(OH)22. Due to no prior information on the 
putative binding modes between MMP-2 and Gd@C82(OH)22, these 
nanoparticles were initially set aside very far away from the protein 
MMP-2 with a distance larger than the typical cutoff distance of van 
der Waals interactions. Following a similar protocol as in our pre- 
vious study of MMP-9"', four Gd@Cs2(OH)22 were again introduced 
at the tetrahedral (Td) corners of the simulation water box, with a 
minimum distance from the MMP-2 protein surface of 25.0 A. 
Multiple nanoparticles were used in order to (i) expedite binding site 
search on MMP-2, (ii) enhance sampling with various initial config- 
urations between MMP-2 and Gd@C82(OH)22, and (iii) simulate the 
potential clustered binding of multiple Gd@C82(OH)22 nanoparticles 
on the protein. As a control, we also examined the binding dynamics 
of the normal fullerenol C82(OH)22 (See details in Methods) . The force 
field parameterization for Gd@C82(OH)22 was adopted from previous 
studies"*" which was developed based on a similar approach 
widely used in standard force field developments. These parameters 
are compatible with data obtained from experiments and high level 
quantum mechanics calculations (see SI text and Fig. S9 of ref 18). 

Overall, both fuUerene derivatives Gd@C82(OH)22 and C82(OH)22 
frequently interact with various sites of MMP-2 while searching for 
the most thermodynamically stable sites (see more details below). 
The nanoparticles were also sometimes observed to cluster together, 
before or after binding onto MMP-2, as occurred in the previous 
simulations with MMP-9"'. This clustering was also observed in 
the atomic force microscopy ( AFM) and synchrotron-radiation small 
angle X-ray scattering (SR-SAXS) experiments, where Gd@C82(OH)22 
can exist as an aggregate with a radius as large as —22 nm'*. This 
indicates that despite the Gd'* ion charge and multiple surface hydro- 
xylation, Gd@C82(OH)22 is still largely hydrophobic. 

MMP-2 remains intact with Gd@C82(OH)22 interaction. The global 
structural change of MMP-2 upon interaction with Gd@C82(OH)22 
and C82(OH)22 was first investigated here. Figure 2 shows the root 
mean square deviation (RMSD) and fluctuation (RMSF) of MMP-2 
with Gd(a)C82(OH)22 and C82(OH)22, respectively. MMP-2 appeared 
quite stable over the simulation time regardless of the interacting ful- 
lerene derivatives. The RMSD grew no more than ~3A in any tra- 
jectories, fluctuating between lA to 3A from the referenced X-ray 
crystal structure, which indicates that the native fold of MMP-2 
remains stable, with its secondary structures of 5 P-strands and 3 a- 
helices all intact, even after frequent contacts with the nanoparticles. 
The similar interaction of Gd@C82(OH)22 was also found in MMP-9, 
where RMSD was stable and smaller than 3 A over —500 ns in multiple 
trajectories'". Similar to the MMP-9 case, we thus excluded the struc- 
tural deformation as a potential inhibitory route for Gd@C82(OH)22 
on MMP-2. This is in contrast to our previous studies with the single 
wall carbon nanotubes (SWCNTs) or graphene nanosheets, where pro- 
tein structures can be significantly distorted/ destioyed by the stiong 
hydrophobic and/or n-n stacking interactions between the SWCNT/ 
graphene and hydrophobic core residues of the proteins'""""'. 

Even though without any noticeable changes in the overall fold of 
MMP-2, the RMSF did indicate that these nanoparticles, particularly 
Gd@C82(OH)22, can sensitively "detect" the surface exposed residues 
in more sophisticated manner (Fig. 2B). The loop regions seem to be 
much more "accommodative" to the nanoparticles, in addition to 
their intrinsic flexibilities. For instance, loops L34 and L45 in the 
upper rim, as well as loops SC and SI' in the lower rim (which are 
involved in the ligand recognition of MMP-2), all display higher 
flexibility to accommodate the binding of these fuUerene derivatives. 
Both Gd@C82(OH)22 and C82(OH)22 were found to have specific 
interactions with MMP-2 at certain sites on the surface, such as near 
the loop L34, where Gd@C82(OH)22 also displayed more selectivity 
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MMP2 195 LFLVAAHEFCHAMGL EHSQDPCALMAP I YTYTKNFRLSQDD I KG I QETYCAS P 250 
MMP9 395 LFLVAAHEFGHALGLDHSSVPEALMYPMYRFTECPPLHKDDVNC I RHLY 443 

Met-turn SV hC 

Figure 1 | Molecular system (MMP-2 and Gd(a)C82(OH)22) and sequence alignment between MMP-2 and MMP-9. (a) Catalytic domain of MMP-2: two 
Zn^* and three Ca^* are depicted in gray and cyan, respectively, and Met-turn and the ligand specificity loop SI' is highlighted with yellow and 
orange, respectively, (b) MMP-2 is superimposed on the catalytic domain of MMP-9 marked in blue with RMSD of 0.97 A. (c) Structure of 
metallofuUerenol Gd@C82(OH)22, where Gd'* is colored in pink, (d) Sequence alignment between MMP-2 and MMP-9 with 64.3% sequence identity. 



hB 



than C82(OH)22 (more below). Interestingly, MMP-9 also showed 
similar patterns with residues at the loop L34 being more reactive to 
Gd@C82(OH)22 than C82(OH)22, possibly due to the fact that they 
share the same fold'". Again, this proposes that Gd@C82(OH)22 
inhibition might be through specific bindings on MMP-2 instead 
of direct structural deformation due to the lack of extreme hydro- 
phobicity as in SWCNT or graphene. 

MMP-2 has more diversified interaction with Gd@C82(OH)22 
than MMP-9. Next, the nanoparticle binding was characterized in 
greater detail by residue-specific contact analysis (Fig. 3). The residue 
contacts were obtained by calculating the fraction of number of contact 
events over the trajectories. A contact event is coimted whenever any 
heavy-atom pair from a target residue and the nanoparticle is within 
5A. Compared with MMP-9, MMP-2 shows a more diverse interaction 
with both Gd@C82(OH)22 and C82(OH)22 at various surface sites 
(cf. Fig. 3 and Fig S2). Overall, 117 and 105 residues of MMP-2 are 
involved in contacts with Gd@C82(OH)22 and C82(OH)22, respectively, 
for at least 1 % of the entire simulation time. By contrast, only 47 and 79 
residues of MMP-9 participated in contact with Gd@C82(OH)22 and 
C82(OH)22, respectively, for the same time duration. 

In addition, both Gd@C82(OH)22 and C82(OH)22 interact with 
much wider range of amino acid types in MMP-2 than in MMP-9. 
Both fuUerene derivatives are effectively interacting with about 17 
different amino acid types in MMP-2, whereas only 1 1 and 14 dif- 
ferent types of amino acids were found in MMP-9 for the interactions 
with Gd@C82(OH)22 and C82(OH)22, respectively (Fig. S3 and Fig. 
S6 of ref. 18). The effective number N^^of amino acid types were 
assessed by N^^ = exp{S), where S is defined by the sum of -p^ In p„ 
over the contact probabilities p^ of 20 amino acids. This clearly 
indicates that MMP-2 and MMP-9 have different interactions with 



Gd@C82(OH)22 (and C82(OH)22)j and Gd@C82(OH)22 is sensitive 
enough to distinguish the two gelatinases even with their highly 
conserved 3-D folds. 

Despite the overall similarity in contact number statistics, detailed 
contact analyses also reveal that there is significant difference between 
Gd@C82(OH)22 and C82(OH)22 in their characteristic bindings onto 
MMP-2. C82(OH)22 seems to have its contacts with MMP-2 relatively 
evenly distributed over the surface, while Gd@C82(OH)22 displays a 
more concentrated, but multiple, hotspots in binding with MMP-2, 
manifested in the loop regions L23, L34 and L45 in the upper rim and 
the ligand specificity SI' loop regions (Fig. 3a). Meanwhile, Gd@ 
C82(OH)22 is also found to directly contact with the catalytic site of 
MMP-2, the highly conserved zinc-binding motif ^"'HExxHxxGxxH""' 
located between the C-termtnus of helix hB and Met-turn. This direct 
contact, however, is significantly reduced for C82(OH)22, which 
implies that the Gd encapsulation is crucial for the nanoparticle to 
have drug efficacy to MMP-2. 

This concentrated multi-hotspot binding of Gd@C82(OH)22 on 
MMP-2 is in great contrast with MMP-9, where a rather dominant 
single binding mode was observed (Fig. S2b)"*. Our previous study 
showed that Gd@C82(OH)22 mostly contacted with the ligand spe- 
cificity SI ' loop of MMP-9, and more importantly, the zinc-catalytic 
site remained intact from direct binding of either Gd@C82(OH)22 or 
C82(OH)22 even in the extended simulations of up to —500 ns. This 
suggested a potential anti-angiogenic mechanism of Gd@C82(OH)22 
through MMP-9, where Gd@C82(OH)22 inhibits the proteolytic activ- 
ity of MMP-9 via aUostericaUy interfering incoming ligands at the 
ligand specificity SI' loop, rather than direcdy blocking the catalytic 
site. However, the multiple hotspots (binding modes) currently found 
for MMP-2 indicate that Gd@C82(OH)22 might inhibit MMP-2 via a 
combined approach with multiple pathways (see more details below). 
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Figure 2 | RMSD and RMSF of MMP-2. (a) RMSD of MMP-2 fluctuates 
between 1 and 3 A for both Gd@C82(OH)22 (up) and C82(OH)22 (down), 
where each trajectories are represented with different colors, (b) RMSF of 
MMP-2 with Gd@C82(OH)22 (up) and C82(OH)22 (down). Five P-sheets 
and three a-helices seem stable, while loops from both upper and lower 
rims are more sensitive to fuUerene derivatives. 



This contact diversity in MMP-2 may be attributed to the sequence 
dissimilarity from MMP-9 especially on the surface. Although 
MMP-2 and MMP-9 share a relatively large portion of amino acids 
from one-dimensional sequence alignment (i.e., 64.3%), they seem to 
have quite disproportionate sequence similarity depending on res- 
idue locations. Indeed, the buried and exposed sites show 79.0% (of 
100 residues) and 38.6% (of 57 residues) sequence similarity, respect- 
ively (here a residue is counted as "exposed" if its solvent- accessible 
surface area is larger than 40% of the total area; otherwise, "buried") 
(Fig. Id). Since the Zn^^-catalylitic domain is largely determined by 
the buried core with many common residues, it is not surprising to 
see both gelatinases MMPs share a very similar 3-D scaffold in their 
catalytic domain (RMSD = 0.97 A). On the other hand, the ligand 
specificity between MMP-2 and MMP-9 is mainly regulated by sur- 
face residues with relatively lower sequence similarity. This surface 
dissimilarity influences their specific interactions with Gd@ 
Cg2(OH)22, resulting in different binding modes and different inhib- 
itory mechanisms. 

Interestingly, previous studies have also shown that the binding 
specificity of Gd@C82(OH)22 is closely related to the surface electro- 
statics of the target proteins For example, MMP-9 creates a 
strong negative field around the Zn^^catalytic site, thus presumably 
facilitating positively charged ligands (Fig. 84). But since the fuU- 
erenol cage of Gd@Cs2(OH)22 is negatively polarized by the Gd^*- 
ion encapsulation, Gd@Cs2(OH)22 is likely prohibited from directly 
contacting the zinc-reactive site. Instead, Gd@C82(OH)22 was driven 
toward the SI' loop, which has numerous neutral or polar residues, 
through specific hydrogen bondings and hydrophobic interaction. 
On the contrary, MMP-2 forms very different electrostatics, with its 



active site surrounded by rather neutral and/or positive electrostatic 
potential (Fig. S4). In addition, the SI' loop region of MMP-2 seems 
more positively polarized by K230, R233 and K190 residues, as com- 
pared to the mostly neutral SI ' loop in MMP-9. This illustrates how 
the surface residues can differentiate the binding characteristics for 
MMP-2 and MMP-9, which in turn influences their native ligand 
recognition as well as their interaction with nanoparticles. Overall, 
these electrostatic potential analyses support that MMP-2 may have 
multiple contacting hotspots with Gd@C82(OH)22, including both 
the Zn^* catalytic site and the ligand specificity SI' loop. 

Gd@C82(OH)22 could inhibit MMP-2 in multiple binding pathways. 

To further illustrate the possible multiple binding modes, we calcu- 
lated the binding free energy (potential of mean force, PMF) for Gd@ 
C82(OH)22 and C82(OH)22 on MMP-2 (Fig. 4). For consistency, we 
adopted the same reaction coordinates (A(-ont> '^min) used in our 
previous study of MMP-9, where A^ont is the contact area between 
the nanoparticle and MMP-2, and d^^^ is the minimum distance 
between the nanoparticle and catalytic Zn^* ion. The binding free 
energy is calculated by WiA^^nt^ d^^^^) = -RT In piA^^^^i, £?„,;„), where 
p(^cont> drain) is the probability of finding the nanoparticle at position 
(^cont> '^min)- As shown in Fig. 4, the binding free energy surface 
clearly discriminates multiple binding modes of Gd(a)C82(OH)22 on 
MMP-2. We identified two specific binding modes, one near the 
catalytic Zn^* ion (i.e., site i in Fig. 4) and one near the SI' loop 
(i.e., site j'i), and other less-specific binding modes (e.g.. Hi, iv, and 
v) relatively distant from the Zn^* ion, depending on the site 
selectivity and potential relevance to the inhibitory function to 
MMP-2. The bmding mode i is located at (285 A^ 6 A) of the 
PMF surface with AG = —4.34 kcal/mol. This corresponds to the 
aforementioned direct binding mode onto the putative ligand binding 
site. The structure (ij) shows that Gd@C82(OH)22 is tightly packed on 
top of the triad of zinc-coordinating histidines (i.e., H201, H205 and 
H211) interacting with the nearby residues (e.g., K89, E210, Y155, 
H166, F168, P90, A169 and L209). As mentioned earlier, this direct 
blocking is possible partly due to the electrostatically favorable local 
environment generated by the zinc-reactive center. Here, Gd@ 
C82(OH)22 seems able to approach closely to the Zn^* ion, which is 
facilitated by the cation-rc interaction'', in addition to the electrostatic 
attraction between the negatively induced fuUerenol cage (i.e., Gd@ 
C82(OH)22[C82(OH)22]'-) and Zn^+ ion*". This mode was absent in 
MMP-9 even with the highly conservative binding site motif (i.e., 
MMP-9: HEFGHALGLDH and MMP-2: HEFGHAMGLEH) due 
to the strong negative electrostatic shielding from nearby acidic 
residues, such as D113, Elll, D182, and D410. Thus, this direct 
binding represents a straightforward inhibitory pathway of Gd@ 
C82(OH)22 to MMP-2. 

Next, Gd@C82(OH)22 was also stabilized near the ligand specifi- 
city SI' loop with AG = —4.89 kcal/mol, as indicated as site ii at 
(205 A"', 10 A) on the PMF surface. Our simulations revealed that 
two representative binding patterns are possible in this energy basin, 
as shown in structures iij and ii2, respectively, which is consistent 
with those in MMP-9 (cf. Fig. S2b). Gd@C82(OH)22 is capable of 
binding in-between the SI' loop with either the SC-loop (i.e., iij) or 
the highly conservative Met-turn (i.e., JJ2). In both structures, Gd@ 
C82(OH)22 can effectively interact with various types of surface 
exposed residues. In mode iij Gd@C82(OH)22 has frequent contacts 
with the residues of both SC (e.g., L164, L185, K190, V191, G192, 
Y193, and S194) and SI' (e.g., 1222, Y223, T227, and Y228) loops, 
potentially interfering the ligand binding at the pave made of these 
two loops leading to the Zn""* -catalytic center. Similarly, Gd@ 
C82(OH)22 in mode ii2 has favorable interactions with residues 
encompassing the Met-turn and SI' loop (e.g., Q213, D214, P215, 
G216, A217, A220, P221, 1222 and Y223), effectively blocking the 
tunnel formed by the Met-turn and SI' loop that also leads to the 
active site. 



SCIENTIFIC REPORTS | 4:4775 | DDI: 1 0. 1 038/srep04775 



4 



A, [ 




100 120 140 160 180 200 22o 240 

Residue index 



180° 




0.2 



CO 



CO 

o 
O 



%8 100 120 



0.15 



0.05 
0 




140 160 180 200 

Residue index 



220 240 



Figure 3 | Residue-specific contacts on MMP-2 with Gd(aC82(OH)22 and Cs2(OH)22- (a) Gd@C82(OH)22 has a favorable contact with MMP-2 
specifically at both the ligand specificity SI ' loop and Zn^* -coordinated active site, (b) C82(OH)22 contacts are reduced in both the SI ' and the active site 
of MMP-2, while interactions with upper rim residues largely remain. 



On the other hand, the energy minima located at a relatively far 
distance (>17 A) from the catalytic zinc ion (i.e., iii, iv, and v) seem 
less specific in terms of the exact locations, partly due to the reaction 
coordinates used here (we adopted the same reaction coordinates of 
MMP-9 for easy comparison, but they might not be ideal). For 
example, although the energy basin iii (235 A^, 18 A) corresponds 
to a quite stable binding free energy AG = —5.08 kcal/mol, it is 
derived from a mixture of multiple different binding modes which 
happen to share same values of these two reaction coordinates - 
contact area and minimum distance from Zn^* ion - thus, the lower 
binding energy could be biased by the choice of the reaction coordi- 
nates as seen in similar protein folding free energy landscapes'". As 
exemplified in the three snapshots selected for each energy minimum 
(i.e., iii, iv, and v), the Gd@C82(OH)22 contacts are highly "diffused" 
over a large area region (with same minimum distance from Zn^^ 
ion) in the upper rim of MMP-2. Thus, even with a seemingly strong 
binding free energy overall, the "diffused" contact areas indicate that 
Gd@C82(OH)22 makes a non-specific interaction with MMP-2 in 
this upper rim, as clearly seen in the contact map in Fig. 3 as well. 
As for the effect of this diffusive large area contact on the native 
function of MMP-2, it is not immediately clear. But it is possible that 
the upper rim part of the interaction might also provide some 



alternative route for allosteric modulation on the enzyme activity'*^ ''''. 
Therefore, in addition to the aforementioned direct binding mode (to 
the Zn^*- catalytic site) with MMP-2, this diffusive mode is another 
contrasting point for the binding patterns of Gd@C82(OH)22 on 
MMP-2 versus MMP-9, where the binding is dominated by the SI' 
loop site. 

These findings also indicate that Gd@C82(OH)22 can have mul- 
tiple "drivers" for its specific bindings onto a target protein"*"''''^^. 
When Gd@C82(OH)22 is far away from the protein without any 
contact, the first landing is guided by the long range electrostatic 
attraction, as shown in many protein and ligand interactions'*". Con- 
sidering the negatively induced fuUerenol cage, Gd@C82(OH)22 is 
likely to be attracted onto the exposed positive residues. For example, 
we found a cluster of positively charged residues around each 
detected Gd@C82(OH)22 binding site: R88, K89, and K91 near the 
binding mode i, and K190, K230, and R233 around the binding mode 
a on MMP-2 (Fig. S4). But once it is closer to the protein, the 
interaction becomes more complicated with multiple forces in play, 
including cation-7t interaction, TZ-n stacking, hydrophobic inter- 
action, hydrogen bonding, as well as electrostatic interactions^'''". 
Consequently, Gd@C82(OH)22 can be driven into a thermodynami- 
cally more stable site, which may be different from its initial contact 
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Figure 4 | Potential of mean force and representative binding modes, (a) Binding free energy surface of Gd@C82(OH)22 reveals several specific binding 
modes relevant to inhibition of MMP-2 (i.e., modes i and ii). (b) Binding free energy surface of C82(OH)22 reveals no specific binding modes related with 
the protein inhibition. In representative binding modes, the Zn^*-coordinated active site is depicted with an orange ball (Zn^*) coordinated with 
three histidines, and contacting residues to Gd@C82(OH)22 are colored depending on residue polarity and charge: blue, red, green and white for positively 
charged, negatively charged, polar (non-charged), and non-polar residues, respectively. 



point. It is noteworthy that the proposed binding sites (i.e., modes i 
and ii) are composed of a number of hydrophobic and aromatic 
residues as well as polar and charged ones. The similar binding 
reaction pathway was also observed in Gd@C82(OH)22's binding 
on MMP-9, where Gd@C82(OH)22 was first attracted to MMP-9 
via non-specific long range electrostatic interaction with K433 and 
R440, followed by more specific hydrophobic and hydrogen bonding 
interactions into the ligand specificity SI' loop'". 

Moreover, we noticed a significant difference between Gd@ 
C82(OH)22 and C82(OH)22 when interacting with MMP-2. While 
both specific and non-specific binding modes are possible for Gd@ 
C82(OH)22, C82(OH)22 is more likely to interact with MMP-2 via 
non-specific modes. For example, the binding mode i of C82(OH)22 
displayed an even larger contact area than that of Gd@C82(OH)22, but 
the structure analysis showed that C82(OH)22 was actually located 
among the three loops (i.e., L23, L34 and L45) in the upper rim with 
arelativelylowbindingfreeenergy(i.e., AG = —4.19 kcal/mol). More 
dramatically, the energy basins ii and Hi, with relatively stable binding 
free energies of —4.84 and —5.14 kcal/mol, respectively, were com- 
posed of various different binding modes, implying the binding sites 
are not well specified (i.e. very diffusive) in each energy basin. Com- 
pared with Gd@C82(OH)22, the residue contacts are largely reduced 
near the critical Zn^*- coordinated catalytic site and the ligand spe- 
cificity S 1 ' loop, as shown in Fig. 3. Interestingly, C82(OH)22 was also 
shown to have a non-specific interaction with MMP-9 with signifi- 
cantly reduced inhibitory effect on its proteolytic function (Fig. S2a) 



Thus, our result is consistent with the in vitro experiments that Gd@ 
C82(OH)22 dramatically reduces proteolytic activity of MMP-2 and 
MMP-9, while the normal fullerenol C82(OH)22 (or C6o(OH)2o) has 
much reduced effect on the enzyme activity of the MMPs. (Due to low 
chemical stability of C82(OH)22, a smaller C6o(OH)2o was used as con- 
trol in our previous experiment. However, we selected C82(OH)22 as 
control in our simulation to better capture the effect of the metal ion 
Gd'* while retaining the molecular size to be comparable with Gd@ 
C82(OH)22. The previous study with both C82(OH)22 and C6o(OH)22 as 
control runs shows similar results"*.) 

These findings might have shed light on the potential inhibitory 
mechanism of Gd@C82(OH)22 on MMP-2. Two specific binding 
modes (i.e., i and ii) were identified for Gd@C82(OH)22's inhibition 
of MMP-2: (i) one for direct blocking of the active site of MMP-2, and 
(ii) the other for exocite binding at the ligand specificity ST loop. In 
the former, Gd@C82(OH)22 directly blocks the catalytic site of MMP- 
2, thus the ligand binding is effectively prohibited or seriously inter- 
rupted, resulting in a decrease in the enzyme activity. In the latter, 
Gd@C82(OH)22 inhibits MMP-2 by allosteric modulation at the SI' 
loop. Although it is indirect, the exocite binding can enhance the 
selectivity of Gd@C82(OH)22 considering the importance of the SI' 
loop in the ligand recognition. Despite the dual inhibitory pathways, 
the overall inhibition capability of Gd@C82(OH)22 toward MMP-2 
can be lower than that toward MMP-9, as indicated in our previous in 
vitro assays (Fig. SI)"*. In the MMP-9 case, Gd@C82(OH)22 predo- 
minandy binds to the S 1 ' loop, resulting in a very strong and exclusive 
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binding, while in the MMP-2 case, multiple contact hotspots with 
many non-specific ones (e.g., mode Hi, ivand vin Fig. 4a) co-exist, 
which might spread the strength thin on each one. 

In MMP inhibitor (MMPi) designs, the zinc-active site has been 
intensively targeted owing to its obvious role in proteolysis. Small 
organic molecules are often designed as potential MMPi's with high 
binding affinities, but almost all of them fail in ultimate FDA 
approval due to the lack of target selectivity and thus resulting in 
serious side effects'". The lack of target selectivity in these MMPi's is 
mainly caused by the highly homologous binding site structures 
among the MMPs''^'''^. Meanwhile, aUosteric modulation was dis- 
cussed as an alternative route to obtain a better target selectivity^'. 
For example, prymidine dicarboxamide derivatives were shown to 
selectively inhibit MMP-13 by binding to the ligand specificity SI' 
loop with no direct effect on the active site'"'. The ligand specificity 
SI' loop is found to be a region with the least similarity in terms of 
length and/or sequence among the available MMPs, thus implying a 
potential target site. Our current simulations show that Gd@ 
C82(OH)22 can potentially inhibit MMP-2 allosterically by interfer- 
ing an incoming ligand at the ligand specificity SI ' loop, in addition 
to its direct binding to the zinc catalytic site. 

Taken together, the allosteric modulation through the ligand spe- 
cificity SI' loop could be a useful and specific inhibitory route for 
metallofuUerenol Gd@C82(OH)22 on both MMP-2 and MMP-9. By 
acting on the SI' loop, Gd@C82(OH)22 is shown to effectively inter- 
fere the native ligand from the putative binding. While this was found 
previously as the most probable pathway for MMP-9, our current 
simulations reveal that MMP-2 can be inhibited by Gd@C82(OH)22 
through both direct binding at the Zn^*- catalytic site and indirect 
binding at the SI' loop site. 

Conclusion 

In this study, we used molecular dynamics simulations to study the 
interaction between the metallofuUerenol Gd@C82(OH)22 and 
matrix metalloproteinases MMP-2 with atomic detail. Although 
Gd@C82(OH)22 was originally developed as a high contrasting agent 
for MRI, it has shown eminent efficacy for various cancers including 
hepatoma, breast cancer, and pancreatic cancer. Among the possible 
anti-tumor mechanisms, Gd@C82(OH)22 was found to inhibit the 
proteolytic activity of MMP-2 and MMP-9, which is critical in down- 
regulating the angiogenesis and tumor extra cellular matrix degra- 
dation, and hence restricting tumor proliferation and metastasis. Our 
results show that Gd@C82(OH)22 can have specific interactions with 
MMP-2, which are closely related to the proteolytic function of 
MMP-2. Similar to MMP-9, Gd@C82(OH)22 also acts on the ligand 
specificity SI' loop, which implies that Gd@C82(OH)22 can induce 
allosteric modulation on MMP-2 proteolysis by residing in-between 
the SI ' loop with either the SC loop or the Met- turn. In addition, Gd@ 
C82(OH)22 can also interact with MMP-2 directly at the zinc-reactive 
site, which was not found in its binding with MMP-9. Thus, Gd@ 
C82(OH)22 achieves its inhibition to MMP-2 via dual binding modes, 
with both direct and indirect bindings. Detailed analyses show that the 
local electrostatic environment differentiates the Gd@Cs2(OH)22 
binding between MMP-2 and MMP-9. The residue constituents on 
each protein surface are quite different, which results in a favorable 
electrostatic environment for Gd@C82(OH)22 at the Zn""* -coordi- 
nated reaction site of MMP-2, but not MMP-9, even though both 
share a close fold with a fairly high overall sequence similarity. 

Various types of nanomedicines have been actively developed in 
recent years in order to complement and/or overcome the deficits of 
traditional small molecule based drugs. Among these efforts, a con- 
siderable amount of researches have been devoted to the develop- 
ment of novel nanocarriers for conventional drug delivery. Our 
current approach focuses on another important route of using nano- 
particles directly as nanodrugs for terminal illness such as cancer. This 
study provides new insights into the de novo design of nanodrugs 



in addition to a better understanding of the various inhibitory 
mechanisms. 

Methods 

In our simulation, we used a 160-residue catalytic domain of human MMP-2 deter- 
mined from X-ray crystallography''''. Four Gd@C82(OH)22 were initially located at the 
tetrahedral corners of the cubic simulation box at least 25 A away from MMP-2, since 
we had no prior information on the putative binding modes of Gd@C82(OH)22 on 
MMP-2. Multiple copies of metaUofuUerenols were employed in order to not only 
facilitate the binding-site searching on MMP-2, but also understand the clustered 
behavior of the nanomolecules in the aqueous media. The composite system was then 
immersed in a cubic 90 A X 90 A X 90 A water box, and neutralized with counter 
ions. Then, the system was ionized with 1 00 mM NaCl, resulting in the total system of 
—69,000 atoms with 22,000 explicit water molecules. Similarly, we configured a 
separate system for Cs2(OH)22 as the control run. The CHARMM22 (c32bl para- 
meter set) force field'^^ and TIP3P water model*^ were used for the protein and water 
molecules, respectively. The nonbonding parameters for Gd@C82{OH)22 and 
C82(OH)22 were separately prepared by the DFT-level quantum mechanical calcu- 
lations, and examined according to the protocol for CHARMM force field develop- 
ment'^. The non-bonding interactions were treated with a typical cutoff of 12 A, while 
the long-range electrostatic interaction were enumerated with the particle-mesh 
Ewald method*". Before production runs, the system was prepared with energy 
minimization for 20,000 steps, followed by 250 ps of equilibration with a 0.5 fs time 
step. For each system, five trajectories were independently generated for at least 
200 ns with a 2 fs time step, where the initial structure was differently configured for 
each run by randomly rotating MMP-2 in order to avoid statistical bias as much as 
possible. AU simulations were carried out at 1 atm and 310 K with the NAMD2 
software massively parallelized for IBM Blue Gene machine^". 
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